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Abstract 



We present a preliminary measurement of the electron angular distribution 
parameter a2 in — > ev events using data collected by the D0 detector 
during the 1994-1995 Tevatron run. We compare our results with next-to- 
leading order perturbative QCD, which predicts an angular distribution of 
(libai cos9*+a2 cos'^6*), where 9* is the angle between the charged lepton and 
the antiproton in the Collins-Soper frame. In the presence of QCD corrections, 
the parameters ai and 02 become functions of p^, the W boson transverse 
momentum. We present the first measurement of 0:2 as a function of . This 
measurement is of importance, because it provides a test of next-to-leading 
order QCD corrections which are a non-negligible contribution to the W mass 
measurement. 
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I. INTRODUCTION 



At the Fermilab Tevatron, which operates at a center of mass energy of ^/i =1.8 TeV, 
W and Z bosons are produced in high energy pp coUisions. In addition to probing elec- 
troweak physics, the study of the production of W and Z bosons provides an avenue to 
explore QCD, the theory of strong interactions. The benefits of using intermediate vector 
bosons to study perturbative QCD are large momentum transfer, distinctive event signa- 
tures, low backgrounds, and a well understood electroweak vertex. In this paper we present 
the measurement of the angular distribution of electrons from W boson decays based on 
data taken by the D0 collider detector during the 1994-1995 Tevatron running period. 

In the parton model, W and Z intermediate vector bosons are produced, at lowest 
order, in head-on collisions of qq constituents of the proton and antiproton, and cannot 
have any transverse momentum. This purely electroweak tree-level process is determined 
by the {V — A) character of electroweak interactions and leads to a rather simple angular 
dependence of the cross section (measured by UAl [Ij-^]): 

"^'^ oc(l±cosr)2 (1) 



d{cos6 

where 6* is the lepton angle in the Collins-Soper rest frame 0] of the W boson in which the 
z-axis is defined as the bisector of the proton momentum and the negative of the antiproton 
momentum. 

At large transverse momentum {pt > 20 GeV), the cross section is dominated by the 
radiation of a single parton with large transverse momentum. Perturbative QCD is 
therefore expected to be reliable in this regime. The additional gluon or quark jet in high 
Pt collisions alters the helicity-state of the W boson; a calculation to next-to-leading order 
in QCD leads to nine helicity amplitudes for the various contributing processes (see 0). 
The cross section can be written in terms of the azimuthal and polar angles. Integrating 
over the azimuthal angle leads to: 

C{l + P{W)aicos9* + a2Cos^9*) (2) 



dqj^ dy d cos 9* 



where P{W) is the polarization of the W boson. The angular parameters ai and a2 are 
functions of the transverse momentum of the W boson, p^ , as shown in figure |I]. In this 
paper we present a measurement of the parameter a2 as a function of pj^. The angular 
parameter ai could not be measured by the D0 detector during Tevatron Run 1 because the 
lack of a magnetic field, for sign identification of electrons, makes it impossible to determine 
the polarization of the W boson. For Tevatron Run II, a central solenoid magnet will be 
installed, allowing for the measurement of both angular parameters. 

The measurement of a2 can serve as a probe for perturbative QCD independent of inclu- 
sive measurements. Since the transverse mass of the W boson is correlated with the decay 
angle of the lepton, QCD effects introduce a systematic shift to the W mass measurement at 
D0 , where a fit to the transverse mass distribution is used to determine the W mass. The 
shift, introduced by perturbative QCD, is 0(40 MeV) for events with px < 15 GeV used 
in the mass measurement. In Run II, when the total error of the W mass will be reduced 
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from the current 105 MeV 0,|TO| to an estimated 50 MeV for 1 fb ^ and to about 30 MeV 
for 10 fb"^ iTTl , a good understanding of this systematic shift is important. 




20 40 60 80 100 120 140 160180 20 40 60 80 100 120 140 160 180 



Pj(W) in GeV p^(W} in GeV 



FIG. 1. The angular parameters 02 (left) and ai (right) as a function of p^^. 



II. EXPERIMENTAL METHOD 
A. Extraction of the lepton angle 

To directly measure the decay angle 6* of the electron in the Collins-Soper frame, all 
momenta in the lab frame have to be known to perform the boost to this particular rest 
frame of the W Boson. This is not possible, however, since the longitudinal momentum 
of the neutrino cannot be measured. A solution to this problem is to use the correlation 
between cos 6* and the transverse W mass to infer cos 6* on a statistical basis. This is done 
using a Bayesian approach. 

Figure H generated from Monte Carlo events that were run through a parameterized 
detector simulation (see [p!^-p^) shows the correlation of the smeared W transverse 
mass and the true value of cos 6*. By correlating the smeared transverse mass as it would 
be measured in the D0 detector and the true (unsmeared) value of cos 6'*, the Bayesian 
analysis described below will yield the unsmeared angular distribution. Note, however, that 
the angular distribution obtained this way is the distribution for accepted events which is 
different from the l+P{W)ai cos 6*+a2 cos^ 6* distribution expected for all events. Since the 
correlation between the angle and the transverse mass depends on the transverse momentum 
of the W boson, a separate correlation plot is used for each px bin. The correlation does 
not depend on the angular parameters ai, and 02, however, as can be seen from equation 
10: 



= — — X 



V2 



2^00 + ^17^ + 027^ - 2(- sin^ 9* + 7^(1 - cos^ 0* sin^ 9*) (3) 
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where is the invariant mass of the eu system and the various parameters are defined as: 





w 
Pt 


7 = 




rrieu 


Oo = 


sin^ 9* 


Oi = 


2 sin^ e* 


02 = 


(1 — cos 



(sin^ 0* — cos' 

2 I* „„„2 n*\2 



,2 J * 



cos'^ 



2 /)*\ 



This equation provides the analytical expression for the dependence of the transverse mass 
on the angles 6* and (p*, where (p* is the azimuthal angle in the CoUins-Soper frame with 



respect to the x-axis: |15]. Note that in this analysis we integrate over 
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COS 



FIG. 2. Smeared W transverse mass versus true cos 9* for pT < 10 GeV. Acceptance cuts are 
applied to this plot. This correlation plot is used to infer the cos 9* distribution from the measured 



rrij!^ distribution. 



To obtain an angular distribution from a measured transverse mass distribution, the trans- 
verse mass distribution has to be inverted using the following probability function: 



/(cos6'*|mT) 



g{mT\ cos^*)/i(cos 9*) 
J g{'mT\cos6*)h{cos6*)dcos6* 



(4) 



where: 



g{mT\ cos6*) is the probability of measuring rriT given a certain cos 6'* value (obtained 
from Monte Carlo) ; 

h{cos 6*) is the prior probability for cos 6*: (1 + cos^ 6*). It refiects all prior knowledge, 
i.e. the angular distribution in the absence of QCD effects; 

/(cos^*|mT) is the probability that an event with transverse mass has a decay 
angle cos 6**. 
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It is important to note that g{mT \ cos e*) = ''J^i'----p^-' where i is the bin number m rriT 
and j is the bin number in cos0*. The resulting cos 6'* distribution is then the angular 
distribution for accepted events. 

The angular distribution can now be inferred from the measured distribution by 
integrating /(cos6'*|mr) over m^: 

all rriT bins 

iV, = Ar™"/(cos^;|mT.) (5) 

i 

where Nj is the number of events in cos^^* bin j. 

To measure the angular parameter 02 a series of cos 9* templates is generated from Monte 
Carlo in bins of , each normalized to unity. Each of the templates is generated taking the 
transverse mass distribution from a high statistics Monte Carlo sample for a specific 02 value 
and converting it into an angular distribution by means of the Bayesian method described 
previously. Detector effects are included by applying smearing, efficiency, and acceptance 
corrections in the Monte Carlo program. The angular distribution obtained from data is 
then compared to these templates to determine which value for a2 fits best. Figure |^ shows 
a series of such templates for < 10 GeV. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cos 

FIG. 3. Templates of the angular distribution for various 02 values for < 10 GeV. These 
templates are obtained from Monte Carlo after acceptance cuts have been applied which results in 
the drop-off at small angles. 



B. Data Selection 

The W sample has been selected from data taken during the 1994-95 run of the Tevatron, 
and corresponds to an integrated luminosity of 84.5 pb~^. 

The measurements of the W boson angular decay distribution used the decay mode 
W — i> eu. Electrons were detected in hermetic, uranium/liquid- argon calorimeters with an 
energy resolution of about 15%/ ^^/E[GeV). The calorimeters have a granularity of Arj x 
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= 0.1 X 0.1, where t] is the pseudorapidity and is the azimuthal angle. Electrons 
were accepted in the region \t]\ < 1.1 (central calorimeter). We assume that the transverse 
momentum of the neutrino is given by the calorimetric measurement of the missing transverse 
energy {^t) in the event. 

Electrons from W and Z boson decays tend to be isolated. Thus, we made the cut 

EtotjOA) - Eem{0.2) 
Eem{0.2) 

where £'tot(0.4) is the energy within AR < 0.4 of the cluster centroid {AR = a/ Arj"^ + A0^) 
and i?£;M(0.2) is the energy in the EM calorimeter within AR < 0.2. 

At trigger level, a single electron with Et greater than 20 GeV was required. Off-line, a 
tighter requirement on the electron quality was introduced to reduce the background level 



from QCD dijet events, especially at high transverse momentum ||T6[. Electron identification 
was based on a likelihood technique. Candidates were first identified by finding isolated 
clusters of energy in the EM calorimeter with a matching track in the central detector. 
We then cut on a likelihood constructed from the following four variables: the from a 
covariance matrix which measures the consistency of the calorimeter cluster shape with that 
of an electron shower; the electromagnetic energy fraction (defined as the ratio of the energy 
of the cluster found in the EM calorimeter to its total energy) ; a measure of the consistency 
between the track position and the cluster centroid; and the ionization dE/dx along the 
track. We require one central isolated electron and > 25 GeV. The event is rejected 
if there is a second electron and the dielectron invariant mass lies in the range 75 — 105 
GeV. A total of 41173 central W candidates passed these cuts. A parametric Monte Carlo 
program was used to simulate the D0 detector response and calculate the kinematic 
and geometric acceptance as a function of ttit and pt- The detector resolutions used in the 
Monte Carlo were determined from data, and were parameterized as a function of energy 
and angle. The relative response of the hadronic and EM calorimeters was studied using the 
transverse momentum of the Z boson as measured by the pt of the two electrons compared 
to the hadronic recoil system in the Z event. 



C. Backgrounds 

In order to use the transverse mass distribution of W Bosons to measure the angular 
distribution, the size of the backgrounds and their dependence on the two variables of interest 
here, transverse momentum and transverse mass, have to be estimated. The following 
sections describe how the four dominant backgrounds are calculated and how they depend 
on transverse mass and momentum. 



1. QCD 

QCD dijet events in which a jet is misidentified as an electron and the energy in the 
event is mismeasured which results in large missing transverse energy pose the largest overall 
background. The reason for this is the very large dijet cross section compared to the W 
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cross section. This background is estimated from a data sample that is dominated by fake 
electrons which can be obtained by replacing the cut for W events with a complementary 
cut $rp < 15 GeV which selects events with fake electrons from multijet events with only a 
small contamination from real W events. The overall QCD background fraction is fn^cjy = 
(0.95±0.6)% {f^cD = (0.77±0.6)%) without (with) a transverse mass cut of 50 < mf < 90 
GeV imposed, which is the range used in the Bayesian analysis. 



2. Z^ee 

Z events can look like W events if one electron is lost in an uninstrumented or under- 
instrumented region of the detector and the resulting energy imbalance amounts to large 

. This background can only be estimated using Monte Carlo Z events. The number of 
such Z events present in the W sample is calculated by applying the W selection cuts to 
HERWIG Z ^ ee events that were passed through a GEANT based simulation of the D0 



detector and were overlaid with events from random pp crossings. The overall background 
fraction is/^ = 0.0091 ± 0.0013 for all pr ■ 



3. Top-Antitop 



The top background is calculated in a similar way as the Z background from HERWIG tt 
events. The overall background fraction isf}^ = 0.0016 ± 0.0005 for all pt ■ 



4. W ^TV 

W ^ T events where the r subsequently decays into an electron and two neutrinos are 
indistinguishable from W ^ ev events. This background is included in the parameterized 
Monte Carlo described above: 2.22 % of the events are generated as r's so that the plot 
correlating cos^* and simply takes these events into account. 

Figure ^ shows the transverse mass distribution of W —>■ eu candidates together with 
all backgrounds (excluding r) in four pjl^ bins. All background rates are rather small. The 
background shapes shown here are fitted with heuristic functions and the errors on the fits 
are used to determine the systematic errors in the measurement of 0:2 due to background 
subtraction. 



D. The measurement of 02 



To obtain the angular distribution for W events from data, the transverse mass distribu- 
tion is inverted according to the Bayesian method described in Section II A . Figure ^ shows 
the background subtracted transverse mass distributions for the four pr bins used in this 
analysis. The x^-values per degree of freedom are 1.67, 1.77, 1.02, and 1.10, respectively, in 
order of increasing pt ■ 
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m^[GeV] m^[GeV] 
FIG. 4. Summary of all backgrounds in four px bins. 



To extract the angular parameter 02 from the angular distribution obtained by inverting 
the transverse mass distribution a log-likelihood method is used, 



all cos 8* bins 

ln£= ^ riilnpi (6) 

i 

where Pi is the normalized population of a cos 6* bin for one of the Monte Carlo templates 
and rii is the population of the same bin in the angular distribution obtained from data. 
The statistical errors for 0:2 are taken at the points where \nC drops by 0.5 units. 

In figure ^the angular distributions obtained from data are compared to the Monte-Carlo 
templates that fit best. In figure |^ the log-likelihood distributions for 02 are shown in the 
four pt ranges. To estimate the sensitivity of this experiment the of 02 distribution 
is calculated with respect to the next-to-leading order QCD prediction and with respect 
to {V — A) theory in the absence of QCD. The with respect to the QCD prediction 
is 0.9/4 dof which corresponds to a probability of 93%. The with respect to the no- 
QCD prediction is 7.0/4 dof which corresponds to 14%. Using the odds-ratio technique, the 
next-to-leading order QCD calculation is found to be preferred by ~ 2cr over the calculation 
in the absence of QCD. The results of this measurement, including the dominant sources 
of error, are summarized in figure ^ and table |. The systematic errors are estimated by 
varying the relevant parameter (like background shape, overall background rate, or modeling 
parameters) in the Monte Carlo by their errors and rerunning the analysis program resulting 
in a varied angular parameter 02- 
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mT,[GeV] 



80 90 

m^[GeV] 




mT.[GeV] 



mT,[GeV] 



FIG. 5. Background subtracted transverse mass distributions(crosses) in four px bins compared 
to Monte Carlo predictions (lines). (D0 preliminary) 



PT 


< pr < 10 


10 < PT < 20 


20 < PT < 35 


35 < PT < 200 


a2 


1.07 ±0.12 


0.81 ±0.23 


0.52 ±0.35 


0.19 ±0.40 


a2, predicted 


0.99 


0.90 


0.69 


0.23 


mean px 


4.4 


12.6 


25.7 


53.6 


error QCD 


±0.022 


±0.048 


±0.086 


±0.069 


error QCD shape 


±0.035 


±0.004 


±0.026 


±0.008 


error Z 


±0.002 


±0.01 


±0.01 


±0.04 


error Z shape 


±0.00 


±0.01 


±0.01 


±0.02 


error top 


±0.0005 


±0.0002 


±0.001 


±0.008 


error hadronic response 


±0.013 


±0.00 


±0.06 


±0.07 


error hadronic recoil model 


±0.04 


±0.09 


±0.10 


±0.10 



TABLE I. Central values and statistical errors for 02 and systematic errors due to backgrounds 
and the hadronic energy scale and resolution. 

III. CONCLUSIONS 

Using data taken with the D0 detector during the 1994-1995 Tevatron collider run, we 
have presented the measurement of the angular distribution of decay electrons from W boson 
events. A next-to-leading order QCD calculation is preferred by ^ 2a over a calculation 
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FIG. 6. Angular distributions for data (crosses) compared to Monte Carlo templates (lines) 
four different pj- bins. (D0 preliminary) 
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FIG. 7. Log likelihood functions for four different pT bins. 
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order QCD calculation by Mirkes (curve) and calculation in the absence of QCD (horizontal line). 



where no QCD effects are taken into account. 
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